A single cell atlas of the mouse seminal vesicle.

During mammalian reproduction, sperm are delivered to the female reproductive tract bathed in a complex medium known as seminal fluid, which plays key roles in signaling to the female reproductive tract and in nourishing sperm for their onwards journey. Along with minor contributions from the prostate and the epididymis, the majority of seminal fluid is produced by a somewhat understudied organ known as the seminal vesicle. Here, we report the first single-cell RNA-seq atlas of the mouse seminal vesicle, generated using tissues obtained from 23 mice of varying ages, exposed to a range of dietary challenges. We define the transcriptome of the secretory cells in this tissue, identifying a relatively homogeneous population of the epithelial cells which are responsible for producing the majority of seminal fluid. We also define the immune cell populations - including large populations of macrophages, dendritic cells, T cells, and NKT cells - which have the potential to play roles in producing various immune mediators present in seminal plasma. Together, our data provide a resource for understanding the composition of an understudied reproductive tissue with potential implications for paternal control of offspring development and metabolism.

An Early Islet Transcriptional Signature Is Associated With Local Inflammation in Autoimmune Diabetes.

Identifying the early islet cellular processes of autoimmune type 1 diabetes (T1D) in humans is challenging given the absence of symptoms during this period and the inaccessibility of the pancreas for sampling. In this article, we study temporal events in pancreatic islets in LEW.1WR1 rats, in which autoimmune diabetes can be induced with virus infection, by performing transcriptional analysis of islets harvested during the prediabetic period. Single-cell RNA-sequencing and differential expression analyses of islets from prediabetic rats reveal subsets of ß- and α-cells under stress as evidenced by heightened expression, over time, of a transcriptional signature characterized by interferon-stimulated genes, chemokines including Cxcl10, major histocompatibility class I, and genes for the ubiquitin-proteasome system. Mononuclear phagocytes show increased expression of inflammatory markers. RNA-in situ hybridization of rat pancreatic tissue defines the spatial distribution of Cxcl10+ ß- and α-cells and their association with CD8+ T cell infiltration, a hallmark of insulitis and islet destruction. Our studies define early islet transcriptional events during immune cell recruitment to islets and reveal spatial associations between stressed ß- and α-cells and immune cells. Insights into such early processes can assist in the development of therapeutic and prevention strategies for T1D.

Intratracheally administered LNA gapmer antisense oligonucleotides induce robust gene silencing in mouse lung fibroblasts.

The lung is a complex organ with various cell types having distinct roles. Antisense oligonucleotides (ASOs) have been studied in the lung, but it has been challenging to determine their effectiveness in each cell type due to the lack of appropriate analytical methods. We employed three distinct approaches to study silencing efficacy within different cell types. First, we used lineage markers to identify cell types in flow cytometry, and simultaneously measured ASO-induced silencing of cell-surface proteins CD47 or CD98. Second, we applied single-cell RNA sequencing (scRNA-seq) to measure silencing efficacy in distinct cell types; to the best of our knowledge, this is the first time scRNA-seq has been applied to measure the efficacy of oligonucleotide therapeutics. In both approaches, fibroblasts were the most susceptible to locally delivered ASOs, with significant silencing also in endothelial cells. Third, we confirmed that the robust silencing in fibroblasts is broadly applicable by silencing two targets expressed mainly in fibroblasts, Mfap4 and Adam33. Across independent approaches, we demonstrate that intratracheally administered LNA gapmer ASOs robustly induce gene silencing in lung fibroblasts. ASO-induced gene silencing in fibroblasts was durable, lasting 4-8 weeks after a single dose. Thus, lung fibroblasts are well aligned with ASOs as therapeutics.

Single cell transcriptomics reveals dysregulated cellular and molecular networks in a fragile X syndrome model.

Despite advances in understanding the pathophysiology of Fragile X syndrome (FXS), its molecular basis is still poorly understood. Whole brain tissue expression profiles have proved surprisingly uninformative, therefore we applied single cell RNA sequencing to profile an FMRP deficient mouse model with higher resolution. We found that the absence of FMRP results in highly cell type specific gene expression changes that are strongest among specific neuronal types, where FMRP-bound mRNAs were prominently downregulated. Metabolic pathways including translation and respiration are significantly upregulated across most cell types with the notable exception of excitatory neurons. These effects point to a potential difference in the activity of mTOR pathways, and together with other dysregulated pathways, suggest an excitatory-inhibitory imbalance in the Fmr1-knock out cortex that is exacerbated by astrocytes. Our data demonstrate that FMRP loss affects abundance of key cellular communication genes that potentially affect neuronal synapses and provide a resource for interrogating the biological basis of this disorder.

Spatial transcriptomic reconstruction of the mouse olfactory glomerular map suggests principles of odor processing.

The olfactory system's ability to detect and discriminate between the vast array of chemicals present in the environment is critical for an animal's survival. In mammals, the first step of this odor processing is executed by olfactory sensory neurons, which project their axons to a stereotyped location in the olfactory bulb (OB) to form glomeruli. The stereotyped positioning of glomeruli in the OB suggests an importance for this organization in odor perception. However, because the location of only a limited subset of glomeruli has been determined, it has been challenging to determine the relationship between glomerular location and odor discrimination. Using a combination of single-cell RNA sequencing, spatial transcriptomics and machine learning, we have generated a map of most glomerular positions in the mouse OB. These observations significantly extend earlier studies and suggest an overall organizational principle in the OB that may be used by the brain to assist in odor decoding.

Simultaneous profiling of multiple chromatin proteins in the same cells.

Methods derived from CUT&RUN and CUT&Tag enable genome-wide mapping of the localization of proteins on chromatin from as few as one cell. These and other mapping approaches focus on one protein at a time, preventing direct measurements of co-localization of different chromatin proteins in the same cells and requiring prioritization of targets where samples are limiting. Here, we describe multi-CUT&Tag, an adaptation of CUT&Tag that overcomes these hurdles by using antibody-specific barcodes to simultaneously map multiple proteins in the same cells. Highly specific multi-CUT&Tag maps of histone marks and RNA Polymerase II uncovered sites of co-localization in the same cells, active and repressed genes, and candidate cis-regulatory elements. Single-cell multi-CUT&Tag profiling facilitated identification of distinct cell types from a mixed population and characterization of cell-type-specific chromatin architecture. In sum, multi-CUT&Tag increases the information content per cell of epigenomic maps, facilitating direct analysis of the interplay of different chromatin proteins.

scRNA-seq of human vitiligo reveals complex networks of subclinical immune activation and a role for CCR5 in Treg function.

Vitiligo is an autoimmune skin disease characterized by the targeted destruction of melanocytes by T cells. Cytokine signaling between keratinocytes and T cells results in CD8+ T cell infiltration of vitiligo lesions, but the full scope of signals required to coordinate autoimmune responses is not completely understood. We performed single-cell RNA sequencing on affected and unaffected skin from patients with vitiligo, as well as healthy controls, to define the role of each cell type in coordinating autoimmunity during disease progression. We confirmed that type 1 cytokine signaling occupied a central role in disease, but we also found that this pathway was used by regulatory T cells (Tregs) to restrain disease progression in nonlesional skin. We determined that CCL5-CCR5 signaling served as a chemokine circuit between effector CD8+ T cells and Tregs, and mechanistic studies in a mouse model of vitiligo revealed that CCR5 expression on Tregs was required to suppress disease in vivo but not in vitro. CCR5 was not required for Treg recruitment to skin but appeared to facilitate Treg function by properly positioning these cells within the skin. Our data provide critical insights into the pathogenesis of vitiligo and uncover potential opportunities for therapeutic interventions.

Hierarchy of signaling thresholds downstream of the T cell receptor and the Tec kinase ITK.

The strength of peptide:MHC interactions with the T cell receptor (TCR) is correlated with the time to first cell division, the relative scale of the effector cell response, and the graded expression of activation-associated proteins like IRF4. To regulate T cell activation programming, the TCR and the TCR proximal interleukin-2-inducible T cell kinase (ITK) simultaneously trigger many biochemically separate signaling cascades. T cells lacking ITK exhibit selective impairments in effector T cell responses after activation, but under the strongest signaling conditions, ITK activity is dispensable. To gain insight into whether TCR signal strength and ITK activity tune observed graded gene expression through the unequal activation of distinct signaling pathways, we examined Erk1/2 phosphorylation or nuclear factor of activated T cells (NFAT) and nuclear factor (NF)-κB translocation in naïve OT-I CD8+ cell nuclei. We observed the consistent digital activation of NFAT1 and Erk1/2, but NF-κB displayed dynamic, graded activation in response to variation in TCR signal strength, tunable by treatment with an ITK inhibitor. Inhibitor-treated cells showed the dampened induction of AP-1 factors Fos and Fosb, NF-κB response gene transcripts, and survival factor Il2 transcripts. ATAC sequencing analysis also revealed that genomic regions most sensitive to ITK inhibition were enriched for NF-κB and AP-1 motifs. Specific inhibition of NF-κB during peptide stimulation tuned the expression of early gene products like c-Fos. Together, these data indicate a key role for ITK in orchestrating the optimal activation of separate TCR downstream pathways, specifically aiding NF-κB activation. More broadly, we revealed a mechanism by which variations in TCR signal strength can produce patterns of graded gene expression in activated T cells.

Gaining Insight into Vitiligo Genetics through the Lens of a Large Epidemiologic Study.

Several epidemiologic studies and GWASs have implicated genetic factors in the pathogenesis of vitiligo. The report by Kim et al. (2020) describes a prospective cohort study from Korea that has the greatest statistical power to date in addressing the epidemiology of vitiligo inheritance. The authors reported the incidence risk ratios in individuals whose first-degree relatives or spouses are affected, providing clear evidence that both genetic and nongenetic factors influence the pathogenesis of vitiligo.

High-Resolution Mapping of Multiway Enhancer-Promoter Interactions Regulating Pathogen Detection.

Eukaryotic gene expression regulation involves thousands of distal regulatory elements. Understanding the quantitative contribution of individual enhancers to gene expression is critical for assessing the role of disease-associated genetic risk variants. Yet, we lack the ability to accurately link genes with their distal regulatory elements. To address this, we used 3D enhancer-promoter (E-P) associations identified using split-pool recognition of interactions by tag extension (SPRITE) to build a predictive model of gene expression. Our model dramatically outperforms models using genomic proximity and can be used to determine the quantitative impact of enhancer loss on gene expression in different genetic backgrounds. We show that genes that form stable E-P hubs have less cell-to-cell variability in gene expression. Finally, we identified transcription factors that regulate stimulation-dependent E-P interactions. Together, our results provide a framework for understanding quantitative contributions of E-P interactions and associated genetic variants to gene expression.

An atlas of cell types in the mouse epididymis and vas deferens.

Following testicular spermatogenesis, mammalian sperm continue to mature in a long epithelial tube known as the epididymis, which plays key roles in remodeling sperm protein, lipid, and RNA composition. To understand the roles for the epididymis in reproductive biology, we generated a single-cell atlas of the murine epididymis and vas deferens. We recovered key epithelial cell types including principal cells, clear cells, and basal cells, along with associated support cells that include fibroblasts, smooth muscle, macrophages and other immune cells. Moreover, our data illuminate extensive regional specialization of principal cell populations across the length of the epididymis. In addition to region-specific specialization of principal cells, we find evidence for functionally specialized subpopulations of stromal cells, and, most notably, two distinct populations of clear cells. Our dataset extends on existing knowledge of epididymal biology, and provides a wealth of information on potential regulatory and signaling factors that bear future investigation.

DolphinNext: a distributed data processing platform for high throughput genomics.

BACKGROUND: The emergence of high throughput technologies that produce vast amounts of genomic data, such as next-generation sequencing (NGS) is transforming biological research. The dramatic increase in the volume of data, the variety and continuous change of data processing tools, algorithms and databases make analysis the main bottleneck for scientific discovery. The processing of high throughput datasets typically involves many different computational programs, each of which performs a specific step in a pipeline. Given the wide range of applications and organizational infrastructures, there is a great need for highly parallel, flexible, portable, and reproducible data processing frameworks. Several platforms currently exist for the design and execution of complex pipelines. Unfortunately, current platforms lack the necessary combination of parallelism, portability, flexibility and/or reproducibility that are required by the current research environment. To address these shortcomings, workflow frameworks that provide a platform to develop and share portable pipelines have recently arisen. We complement these new platforms by providing a graphical user interface to create, maintain, and execute complex pipelines. Such a platform will simplify robust and reproducible workflow creation for non-technical users as well as provide a robust platform to maintain pipelines for large organizations. RESULTS: To simplify development, maintenance, and execution of complex pipelines we created DolphinNext. DolphinNext facilitates building and deployment of complex pipelines using a modular approach implemented in a graphical interface that relies on the powerful Nextflow workflow framework by providing 1. A drag and drop user interface that visualizes pipelines and allows users to create pipelines without familiarity in underlying programming languages. 2. Modules to execute and monitor pipelines in distributed computing environments such as high-performance clusters and/or cloud 3. Reproducible pipelines with version tracking and stand-alone versions that can be run independently. 4. Modular process design with process revisioning support to increase reusability and pipeline development efficiency. 5. Pipeline sharing with GitHub and automated testing 6. Extensive reports with R-markdown and shiny support for interactive data visualization and analysis. CONCLUSION: DolphinNext is a flexible, intuitive, web-based data processing and analysis platform that enables creating, deploying, sharing, and executing complex Nextflow pipelines with extensive revisioning and interactive reporting to enhance reproducible results.

HIV-1-induced cytokines deplete homeostatic innate lymphoid cells and expand TCF7-dependent memory NK cells.

Human immunodeficiency virus 1 (HIV-1) infection is associated with heightened inflammation and excess risk of cardiovascular disease, cancer and other complications. These pathologies persist despite antiretroviral therapy. In two independent cohorts, we found that innate lymphoid cells (ILCs) were depleted in the blood and gut of people with HIV-1, even with effective antiretroviral therapy. ILC depletion was associated with neutrophil infiltration of the gut lamina propria, type 1 interferon activation, increased microbial translocation and natural killer (NK) cell skewing towards an inflammatory state, with chromatin structure and phenotype typical of WNT transcription factor TCF7-dependent memory T cells. Cytokines that are elevated during acute HIV-1 infection reproduced the ILC and NK cell abnormalities ex vivo. These results show that inflammatory cytokines associated with HIV-1 infection irreversibly disrupt ILCs. This results in loss of gut epithelial integrity, microbial translocation and memory NK cells with heightened inflammatory potential, and explains the chronic inflammation in people with HIV-1.

Recovery of viable endocrine-specific cells and transcriptomes from human pancreatic islet-engrafted mice.

Human pancreatic islets engrafted into immunodeficient mice serve as an important model for in vivo human diabetes studies. Following engraftment, islet function can be monitored in vivo by measuring circulating glucose and human insulin; however, it will be important to recover viable cells for more complex graft analyses. Moreover, RNA analyses of dissected grafts have not distinguished which hormone-specific cell types contribute to gene expression. We developed a method for recovering live cells suitable for fluorescence-activated cell sorting from human islets engrafted in mice. Although yields of recovered islet cells were relatively low, the ratios of bulk-sorted ß, α, and δ cells and their respective hormone-specific RNA-Seq transcriptomes are comparable pretransplant and posttransplant, suggesting that the cellular characteristics of islet grafts posttransplant closely mirror the original donor islets. Single-cell RNA-Seq transcriptome analysis confirms the presence of appropriate ß, α, and δ cell subsets. In addition, ex vivo perifusion of recovered human islet grafts demonstrated glucose-stimulated insulin secretion. Viable cells suitable for patch-clamp analysis were recovered from transplanted human embryonic stem cell-derived ß cells. Together, our functional and hormone-specific transcriptome analyses document the broad applicability of this system for longitudinal examination of human islet cells undergoing developmental/metabolic/pharmacogenetic manipulation in vivo and may facilitate the discovery of treatments for diabetes.

Frequent Loss of IRF2 in Cancers Leads to Immune Evasion through Decreased MHC Class I Antigen Presentation and Increased PD-L1 Expression.

To arise and progress, cancers need to evade immune elimination. Consequently, progressing tumors are often MHC class I (MHC-I) low and express immune inhibitory molecules, such as PD-L1, which allows them to avoid the main antitumor host defense, CD8+ T cells. The molecular mechanisms that led to these alterations were incompletely understood. In this study, we identify loss of the transcription factor IRF2 as a frequent underlying mechanism that leads to a tumor immune evasion phenotype in both humans and mice. We identified IRF2 in a CRISPR-based forward genetic screen for genes that controlled MHC-I Ag presentation in HeLa cells. We then found that many primary human cancers, including lung, colon, breast, prostate, and others, frequently downregulated IRF2. Although IRF2 is generally known as a transcriptional repressor, we found that it was a transcriptional activator of many key components of the MHC-I pathway, including immunoproteasomes, TAP, and ERAP1, whose transcriptional control was previously poorly understood. Upon loss of IRF2, cytosol-to-endoplasmic reticulum peptide transport and N-terminal peptide trimming become rate limiting for Ag presentation. In addition, we found that IRF2 is a repressor of PD-L1. Thus, by downregulating a single nonessential gene, tumors become harder to see (reduced Ag presentation), more inhibitory (increased checkpoint inhibitor), and less susceptible to being killed by CD8+ T cells. Importantly, we found that the loss of Ag presentation caused by IRF2 downregulation could be reversed by IFN-stimulated induction of the transcription factor IRF1. The implication of these findings for tumor progression and immunotherapy are discussed.

Diverse repertoire of human adipocyte subtypes develops from transcriptionally distinct mesenchymal progenitor cells.

Single-cell sequencing technologies have revealed an unexpectedly broad repertoire of cells required to mediate complex functions in multicellular organisms. Despite the multiple roles of adipose tissue in maintaining systemic metabolic homeostasis, adipocytes are thought to be largely homogenous with only 2 major subtypes recognized in humans so far. Here we report the existence and characteristics of 4 distinct human adipocyte subtypes, and of their respective mesenchymal progenitors. The phenotypes of these distinct adipocyte subtypes are differentially associated with key adipose tissue functions, including thermogenesis, lipid storage, and adipokine secretion. The transcriptomic signature of "brite/beige" thermogenic adipocytes reveals mechanisms for iron accumulation and protection from oxidative stress, necessary for mitochondrial biogenesis and respiration upon activation. Importantly, this signature is enriched in human supraclavicular adipose tissue, confirming that these cells comprise thermogenic depots in vivo, and explain previous findings of a rate-limiting role of iron in adipose tissue browning. The mesenchymal progenitors that give rise to beige/brite adipocytes express a unique set of cytokines and transcriptional regulators involved in immune cell modulation of adipose tissue browning. Unexpectedly, we also find adipocyte subtypes specialized for high-level expression of the adipokines adiponectin or leptin, associated with distinct transcription factors previously implicated in adipocyte differentiation. The finding of a broad adipocyte repertoire derived from a distinct set of mesenchymal progenitors, and of the transcriptional regulators that can control their development, provides a framework for understanding human adipose tissue function and role in metabolic disease.

HLA Class II Antigen Processing and Presentation Pathway Components Demonstrated by Transcriptome and Protein Analyses of Islet ß-Cells From Donors With Type 1 Diabetes.

Type 1 diabetes studies consistently generate data showing islet ß-cell dysfunction and T cell-mediated anti-ß-cell-specific autoimmunity. To explore the pathogenesis, we interrogated the ß-cell transcriptomes from donors with and without type 1 diabetes using both bulk-sorted and single ß-cells. Consistent with immunohistological studies, ß-cells from donors with type 1 diabetes displayed increased Class I transcripts and associated mRNA species. These ß-cells also expressed mRNA for Class II and Class II antigen presentation pathway components, but lacked the macrophage marker CD68. Immunohistological study of three independent cohorts of donors with recent-onset type 1 diabetes showed Class II protein and its transcriptional regulator Class II MHC trans-activator protein expressed by a subset of insulin+CD68- ß-cells, specifically found in islets with lymphocytic infiltrates. ß-Cell surface expression of HLA Class II was detected on a portion of CD45-insulin+ ß-cells from donors with type 1 diabetes by immunofluorescence and flow cytometry. Our data demonstrate that pancreatic ß-cells from donors with type 1 diabetes express Class II molecules on selected cells with other key genes in those pathways and inflammation-associated genes. ß-Cell expression of Class II molecules suggests that ß-cells may interact directly with islet-infiltrating CD4+ T cells and may play an immunopathogenic role.

DEBrowser: interactive differential expression analysis and visualization tool for count data.

BACKGROUND: Sequencing data has become a standard measure of diverse cellular activities. For example, gene expression is accurately measured by RNA sequencing (RNA-Seq) libraries, protein-DNA interactions are captured by chromatin immunoprecipitation sequencing (ChIP-Seq), protein-RNA interactions by crosslinking immunoprecipitation sequencing (CLIP-Seq) or RNA immunoprecipitation (RIP-Seq) sequencing, DNA accessibility by assay for transposase-accessible chromatin (ATAC-Seq), DNase or MNase sequencing libraries. The processing of these sequencing techniques involves library-specific approaches. However, in all cases, once the sequencing libraries are processed, the result is a count table specifying the estimated number of reads originating from each genomic locus. Differential analysis to determine which loci have different cellular activity under different conditions starts with the count table and iterates through a cycle of data assessment, preparation and analysis. Such complex analysis often relies on multiple programs and is therefore a challenge for those without programming skills. RESULTS: We developed DEBrowser as an R bioconductor project to interactively visualize every step of the differential analysis, without programming. The application provides a rich and interactive web based graphical user interface built on R's shiny infrastructure. DEBrowser allows users to visualize data with various types of graphs that can be explored further by selecting and re-plotting any desired subset of data. Using the visualization approaches provided, users can determine and correct technical variations such as batch effects and sequencing depth that affect differential analysis. We show DEBrowser's ease of use by reproducing the analysis of two previously published data sets. CONCLUSIONS: DEBrowser is a flexible, intuitive, web-based analysis platform that enables an iterative and interactive analysis of count data without any requirement of programming knowledge.

Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus.

Eukaryotic genomes are packaged into a 3-dimensional structure in the nucleus. Current methods for studying genome-wide structure are based on proximity ligation. However, this approach can fail to detect known structures, such as interactions with nuclear bodies, because these DNA regions can be too far apart to directly ligate. Accordingly, our overall understanding of genome organization remains incomplete. Here, we develop split-pool recognition of interactions by tag extension (SPRITE), a method that enables genome-wide detection of higher-order interactions within the nucleus. Using SPRITE, we recapitulate known structures identified by proximity ligation and identify additional interactions occurring across larger distances, including two hubs of inter-chromosomal interactions that are arranged around the nucleolus and nuclear speckles. We show that a substantial fraction of the genome exhibits preferential organization relative to these nuclear bodies. Our results generate a global model whereby nuclear bodies act as inter-chromosomal hubs that shape the overall packaging of DNA in the nucleus.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon.

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.